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Abstract 

We consider constraints on long-lived charged scalar leptons Z ± in supersymmet- 
ric models, where gravitino is the lightest superparticle. We study the decay and 
hadronization processes of I . We also take into account the significant enhancement 
of 6 Li production due to the formation of the bound-state ( 4 HeZ~); for this purpose, 
we use the reaction rate given by the most recent calculation based on coupled-channel 
method. 



Cosmology provides important constraints on physics beyond the standard model. One of 
the important examples is constraints on long-lived particles from big-bang nucleosynthesis 
(BBN); since the theoretical predictions of the standard BBN scenario are in reasonable 
agreements with observation of light-element abundances, we obtain constraints on long- 
lived particles whose decay may induce photo- and hadro-dissociation processes. 

Here, as a physics beyond the standard model, we consider low energy supersymmetry 
(SUSY), which is strongly motivated as a prominent solution to serious problems in the 
standard model, like the hierarchy problem, the dark-matter problem, and so on. In the 
framework of SUSY models, we assume that one of the charged sleptons is the lightest 
particle among the superpaticle in the minimal-supersymmetric-standard-model (MSSM) 
sector, and that the gravitino is the lightest superparticle (LSP). Such a mass spectrum is 
easily realized; probably the most famous example is gauge-mediated SUSY breaking models 
PP. Thus, such a scenario has been attracted much attention, in particular, in the context of 
gravitino LSP [21EII11EJEJEIEJIH]- Then, the lightest slepton I is usually long-lived and the 
decay of I may spoil the success of BBN. Indeed, the BBN scenario with long-lived charged 
sleptons has been discussed in several works [6J [7J [10] . 

In this letter, we reconsider effects of long-lived charged slepton on BBN. We perform 
a detailed study on the decay processes of slepton and the hadronization processes of the 
decay products. In addition, we also take into account the effect of bound-state formation; 
recently, it has been pointed out that long-lived charged slepton may form bound states with 
light elements [HI [12], [13], which may significantly change the BBN reaction rates [IT] . In 
particular we take a serious attitude towards the bound-state effect by ( 4 HeZ~) pointed out 
by [UJ (which is called /"-catalyzed process) and include it into our analysis, taking into 
account the result of the most recent precise calculation of the bound-state effect [TJ]. 

Let us first summarize the basic properties of the lightest slepton / and the gravitino 
(which is denoted as ip^ hereafter). Since the gravitino is the only superparticle lighter than 
I, the dominant decay mode of I is / — > lip^. The decay rate for this process is given by 



where m z ~ and 777.3/2 are slepton mass and gravitino mass, respectively, and M* ~ 2.4 x 
10 18 GeV is the reduced Planck scale. When the gravitino is much lighter than I, we obtain 



For the case where the next-to-the-lightest superparticle (NLSP) is selectron e, the daughter 
particles are stable. For other cases, on the contrary, the final-state lepton also decays; such 
decay process is sometimes important. In particular, when stau is the NLSP, r-lepton in the 
final state decays subsequently and produces pion(s). As we will discuss later, such pions 
become a source of p <-> n conversion. In addition, when fi- and r-lepton decays, neutrinos 
carry away some amount of energy, which affects the total amount of visible energy emitted 
by the decay of I. 
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In our analysis, we study the decay processes of unstable leptons (in particular, that of 
r-lepton) using PYTHIA package [153 an d derive energy distributions of the decay products. 
For the decay of r-lepton, for example, the averaged number of charged pions produced by 
the decay of r-lepton is 0.94, which changes the 4 He abundance as we will see later. With 
the result given by PYTHIA analysis, we also calculate the (averaged) visible energy emitted 
from I taking into account the energy loss due to the neutrino emission. 

Charged sleptons may also^ decay into 3- and/or 4-body final states as I — > lip^Z^, 



followed by the decay of ^(*) J#jJ Although 3- and 4-body decay processes are sub-dominant, 
they are important for the study of hadro-dissociation processes. When Z^*> decays into qq 
pair, protons, neutrons, and pions are emitted after hadronizing final-state quarks. These 
hadrons induce hadro-dissociation and p <-> n conversion processes, which affect abundances 
of light elements. (For details, see, for example, [IS].) In order to study the effects of hadro- 
dissociation and p ^ n conversion processes, it is necessary to know the hadronic branching 
ratio of I as well as the spectra of final-state hadrons. We have also calculated the decay rate 
for the processes I — > l^^qq. Detailed information about the 3- and 4-body decay processes 
can be found in [10]; we found that our results are consistent with those shown in [10J. We 
have also calculated the invariant-mass and energy distributions of final-state quark pair. 
Then, using PYTHIA package, we calculate energy distributions of p, n, and n^. 

Once the decay rate, total visible energy emitted by the decay of /, and spectra of p, n 
and 7r ± are obtained, we calculate the light-element abundances taking account of photo- 



and hadro-dissociation processes as well as p <-> n conversion processes^! we follow the 
procedure given in [16j in order to calculate the abundances of D, 3 He, 4 He, 6 Li and 7 Li for 
processes except for ( 4 HeZ~) + D — > 6 Li + l~. 

One of the important new ingredients in the current study is the /"-catalyzed process for 



6 Li production [lljcfl In the standard BBN scenario, the production rate of 6 Li is very tiny 
because the El transition process for D( 4 He, 7) 6 Li process is extremely suppressed, resulting 
in very small primordial abundance of 6 Li which is consistent with observation. Recently, 
however, it was claimed that the 6 Li abundance may become much larger than the standard 
BBN prediction if there exists a long-lived heavy charged particle [11]. This is due to the 
enhancedment of the process 

( 4 Her) + D^ 6 Li + r, (3) 

where ( 4 He/~) denotes the bound state formed by 4 He and l~. Importantly, for this pro- 
cess, accurate calculation for the reaction rate is performed with the use of coupled-channel 
method which is known to accurately describe other three-body systems [14J. Thus, we 
adopt the results of [H] and take into account the bound-state effect in the calculation of 
6 Li abundance. 



# 1 Here, we assume that the sleptons are right-handed. In such a case, 3- and 4-body decay processes via 
the emission of are irrelevant. 
# 2 We have omitted contributions from Kaons as discussed in [IB] . 
#3 See also [T^l [T2] for recent progress in this field. 
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In order to estimate the number density of ( 4 He/ ) bound state n^s, we solve the Saha- 



type equation [T7] 



#1 



(rj~i v — 3/2 
2 " G j e EB/T (n4 Uc -n BS )(n h -n BS ), (4) 

where E& ~ 337.33 keV [14] is the binding energy, and n4 He and are total number 
densities of 4 He and l~ (which include those in the bound state), respectively. (In this 
letter, nx is used for the number density of particle X.) Here, we assume that the process 
4 He + ( 4 He/~) is in chemical equilibrium. Then, the total reaction rate R for 6 Li 

production by D- 4 He collision is calculated as 

R=!^R w +(l-!**)R gBBS , (5) 



nine V rn n , 

where Rbs is given by Eq. (4.3) in [Uj and -Rsbbn is the standard-BBN value. 

Once the light-element abundances are theoretically calculated, we compare them with 
observational constraints and derived 95 % C.L. constraints. As observational constraints 
on the primordial abundances of light elements, we adopt those used in [19] except for Y p 
and (n3 H e/n D )p: 

(n D /n H ) P = (2.82 ± 0.26) x 10" 5 , (6) 
(n3 H e/n D )p < 0.83 + 0.27, (7) 
Y p = 0.2516 ±0.0040, (8) 
log 10 (n7 Li /n H ) P = -9.63 ±0.06 ±0.3. (9) 
(m Li /ri7 U ) p < 0.046 ± 0.022 ± 0.084, (10) 

where the subscript "p" is for primordial value (just after BBN), and Y p is the primordial 
mass fraction of 4 He. For the center value of Y p , we have adopted a quite recent value reported 



in [20] in which the authors used new data of Hel emissivitescj and conservatively added 
larger errors of 0.0040 as shown in [22]. For (n3n e /n-D) P we have used most newly-reported 
values of D and 3 He abundances observed in the protosolar cloud [23] . 

Now, we are at the position to show our results. First, we consider the case where 
there is no entropy production after the freeze out of I. In this case, the abundance of I is 
thermally determined; then, for right-handed slepton, the yield variable (which is defined as 
Yj = [nj/s]t<^ T p with s being entropy density) is estimated as [24] 

Fr^7xl0- 14 x (11} 



' " V100 GeV 

Notice that, in the solving Eq. (0J, nj- is given by 



\sYje-^. (12) 
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Figure 1: Constraints on 1113/2 vs. rrif R plane. The dashed line indicates the constraint from 
(n6 Li /n7 Li ) p without taking account of the /"-catalyzed process. We have shaded the region 
where tr becomes lighter than gravitino. 

since ni- = n ; ~+ = 0.5n z ~. 

With the thermal abundance, we perform our analysis for the case where / is tr. The 
result is shown in Fig. HJ For comparison, we also show the result without including the 
/"-catalyzed process (J3j) . As expected, constraints from light elements other than 6 Li is 
unaffected by the inclusion of the /"-catalyzed process ([3]); those constraints are more or less 
consistent with the results obtained in earlier analysis [El [TJ \10\. As one can see, however, 
significant amount of the parameter region is excluded with the overproduction of 6 Li once 
the /"-catalyzed process is considered. For the case where mj = 100 GeV, for example, the 
upper bound on the gravitino mass becomes ~ 0.2 GeV, which was ~ 7 GeV before the 
process ([3]) is taken into account. 

So far, we have assumed that the primordial abundance of / is thermally determined. 
With non-standard evolution of the universe, however, Yj may deviate from the thermal value 
given in Eq. ffTTT) . Thus, we also treat the primordial abundance of / as a free parameter and 
derive upper bound on Yj. For f fi -NLSP, the result is shown in Fig. [21 taking m fR = 100 GeV. 
As one can see, the /"-catalyzed process significantly changes the upper bound; for nif R = 

# 4 For more detailed discussion, see [TJl [18] . 

#5 See also recent value of Y p reported in [21] where the authors adopted larger errors (0.0028) than that 
of [20]. 
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Figure 2: Upper bounds on mf R Y fR as functions of m 3 / 2 for m fR = 100 GeV. The solid 
lines are bounds from various light-element abundances, while the dashed line is the upper 
bound from (neu/n7 Li ) p without taking account of the /"-catalyzed process. Thick dotted 
line shows the contour where the density parameter of the gravitino produced by the decay 
of stau becomes as large as that of cold dark matter. 

100 GeV, for example, the very stringent constraint is obtained when the gravitino mass 
is larger than O(0.1 GeV), which corresponds to r^O(10 3 sec). We have also considered 
the case with m? R = 300 GeV, and the result is shown in Fig. [3J For m? R = 300 GeV, the 
constraints from the hadro-dissociation processes are more stringent than those for m,f R = 
100 GeV. This is mainly because, for the case of mf R = 300 GeV, the decay mode tr — > rZ 
becomes more effective, resulting in larger number of p and n. 

We have also studied another cases that the NLSP is the right-handed selectron cr. In 
this case, the constraints slightly change. First, as we mentioned earlier, in the dominant 
decay mode, tR decays into gravitino and electron, which are both stable. Thus, contrary 
to the f/j-case, the pion production is only via three- and/or four-body decay processes, and 
hence is significantly suppressed. Consequently, the constraint from the p <-> n conversion is 
weakened. Second, the visible energy emitted by the slepton decay is larger for the selectron- 
NLSP case since there is no energy loss due to the neutrino emission (in the two-body decay 
process). As a result, the constraint from the photo-dissociation process becomes more 
stringent than the stau-NLSP case. We checked that the upper bound on Yj from the ratio 
(n3 He /'^D)p becomes smaller by factor ~ 3 compared to the stau-NLSP case. 

Finally, we comment on the bound-state effects on other processes. In our analysis, we 
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Figure 3: Same as Fig. EJ except for m, 



300 GeV. 



consider only the /"-catalyzed process for the 6 Li production. In fact, reaction rates of other 
processes for which the El transition is effective (like T( 4 He,7) 7 Li and 3 He( 4 He, 7) 7 Be) may 
be also enhanced by the bound-state formation. Unfortunately, for those processes, we could 
not find reliable calculation of the reaction rates. However, we expect that the constraint 
from 7 Li is less stringent than that from 6 Li even if the bound-state effect is taken into 
account for Li and 7 Be productions; there are several reasons why we consider so. First, 
observationally, 6 Li is rarer than 7 Li, as one can see in Eq. ( 1101) . Thus, the constraint on 
the model is usually much more sensitive to the change of 6 Li rather than that of 7 Li. In 
addition, reaction rates for the processes T(( 4 He/~), /~) 7 Li and 3 He(( 4 He/~), Z~) 7 Be should 
acquire extra suppressions compared to the process ([3]); for the formar, one should keep in 
mind that the number density of T during BBN is about two orders of mangnitude smaller 
than that of D while, for the latter, the Coulomb suppression becomes much more significant 
compared to the process ([3]). Furthermore, we have also performed numerical calculation by 
using the enhancement factor presented in [11] for 7 Li and 7 Be productions. (Notice that, 
for 6 Li production process ([3]), the formula given in (TTJ predicts reaction rate which is one 
order of magnitude larger than the result of the detailed calculation given in [14J.) Then 
we have checked that, even with these modifications of reaction rates, constraint from 7 Li is 
much less stringent than that from 6 Li. Thus, we expect that the bound-state effects on the 
7 Li and 7 Be production processes do not change our main conclusion. 

Note Added: After the completion of this work, we found the paper [25J which discusses 
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possible effects of the formation of bound states ( 7 LiZ~) and ( 7 Be/~). Such bound states, 
however, do not have significant effect on 6 Li formation and hence our main conclusion is 
unchanged. 
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